Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Mass transfer from the donor of GRS 1915+105 
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Abstract. A scenario for a periodic filling and emptying of the accretion disc of the microquasar GRS1915+105 
is proposed, by estimating the mass transfer rate from the donor and comparing it with the observed accretion 
rate onto the primary black hole. The mass of the Roche-lobe-filling donor (1.2 ±0.2 Mq), the primary black hole 
mass ( 14 ±4 Mq) and the binary orbital period of 33.5 d (Greiner et al. 2001b) predict for the donor spectral 
type and K-magnitude around K6 III and -2.6, respectively. The He-core of 0.28 Mq of such a giant leads to 
evolutionary expansion along the giant branch with a conservative mass transfer rate of Ma = (1-5 ±0.5) x 10~ 8 
Mg/year. On the other hand, the average observed accretion rate onto the primary is ten times larger: M t, s = 2.0 
x (?7/0.1) — 1 (d/12.5kpc) 2 10 -7 Mq/v, where r/ is the efficiency of converting accretion into radiation. We propose 
a duty cycle with (5-10)(?7/0.1) per cent active ON-state. The timescale of the (recurrent) OFF-state is identified 
as the viscosity time scale at the circularization radius (14 Rq) and equals t v i sc = 370 (q/0.001)~ 4//5 years, where 
a is the viscosity parameter in the a-prescription of a classical disc. If the viscosity at the outer edge of the disc 
is small and r/ is close to the maximum available potential energy (per rest mass energy) at the innermost stable 
orbit, the present activity phase may still last another 10 - 20 years. We also discuss other solutions allowing a 
broader range of donor masses (0.6 - 2.4 Mq). 
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1. Introduction 

Greiner et al. (2001a) identified the mass-donating sec- 
ondary star of GRS 1915+105 to be a K-M III giant, in- 
dicating that this prototype microquasar is a low-mass 
' X-ray binary (LMXB). Further, using the Very Large 
Telescope (VLT) and the band-heads of 12 CO and 13 CO , 
Greiner et al. (2001b) managed to obtain the radial veloc- 
ity curve of the secondary. The orbital period of 33.5 days, 
the large mass function f(M) = 9.5 ±3.0 Mq and known 

o 

jet-inclination (70 ) permitted to constrain the primary 
black hole mass between Mbh = (10 - 18) Mq, assuming 
the donor mass to lie between M^ = (1.0 - 1.4) Mq. The 
large BH mass points to rapid rotation since the small- 
est inner disc radii modelled (see e.g. Belloni et al. 1997; 
Vilhu et al. 2001) are as small as 20 km, close to the last 
marginally stable orbit (0.5 R g ) of an extreme prograde 
Kerr- hole of 14 Mq. 

In the present paper we estimate the mass transfer rate 
from the evolving donor but allowing a broader range (0.6 
- 2.4) Mq for its mass to include e.g. a possible stripped 
giant. Further, we estimate the viscosity time scale at the 
circularization radius and the amount of mass accumu- 
lated there. Using the mean observed accretion rate over 
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the past 6 years (via luminosity conversion) we arrive at an 
estimate for the timescale of the possible duty cycle, rele- 
vant also when discussing a link to ultraluminous sources 
(ULX) in other galaxies (King et al. 2001). 



2. Mass transfer rate from the evolving donor 

We assume that the donor fills its Roche lobe and that the 
mass loss is determined by evolutionary expansion along 
the giant branch, conserving the orbital angular momen- 
tum. The properties along the giant branch (luminosity 
and radius) depend mainly on the He-core mass and less 
on the envelope mass. In this case an analytical simplifica- 
tion is possible (Webbink, Rappaport and Savonije, 1983) 
and the procedure is also presented by Verbunt and van 
den Heuvel (1995). In particular, the radius and luminos- 
ity can be fitted with 3rd order polynomials on the core 
mass. 

The growth of the core mass, resulting in an increase 
of the radius, is determined by the luminosity due to hy- 
drogen shell burning which, in turn, depends completely 
on the core mass. In the conservative case, fixing the bi- 
nary parameters and forcing the secondary to fill its Roche 
lobe, it is rather simple to compute the core mass and 
consequently the mass loss from the donor (we use Pop. I 
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Fig. 1. Conservative mass transfer rates from the evolving K-giant donor of GRS 1915+105 using the analytic methods 
by Webbink, Rappaport and Savonije (1983) and Pop. I abundances Z = 0.02. The curves were computed for different 
donor masses (in Mq) using three mass function values inside the range f(M) = 9.5 ±3.0 Mq (Greiner at al. 2001b) 
with P or b fixed to 33.5 d (solid line: f(M) = 9.5, dotted lines: f(M) = 6.5 and 12.5). The dashed line represents the 
analytic expression given by King et al. (2001, their eq. 12). The box shows the mass range suggested by Greiner et 
al. (2001b). 



abundances Z = 0.02; for details see Webbink et al. (1983) 
and Verbunt and van den Heuvel (1995), p. 482). 

The first line in Table 1 gives the results for the best-fit 
masses given by Greiner et al. (2001b) (14 M + 1.2 M , 
P=33.5 d). The second line gives the parameter ranges 
if the donor mass is varied between (0.6 - 2.4) Mq 
and satisfying the mass function constraint f(M)=9. 5+3.0 
Mq. The BH mass varies within these domains between 
9-19Mq. The mass transfer rate depends mainly on the 
donor mass as shown explicitly in Fig. 1. It follows closely 
the analytic expression given by King et al. (2001; their 
eq. 12, the dashed line in Fig. 1). In particular, the un- 
certainties in the mass function have minor effects on this 
relation. 



3. Duty-cycle time scales 

In the conservative case, the mass leaving the donor via 
the Ll-point settles down into a Keplerian orbit around 



the primary BH, the radius of which is called the 'circu- 
larization radius'. This radius is given in Table 1 as com- 
puted from the analytic approximation to numerical data 
(Frank, King and Rainc, 1992 (FKR), p. 56, eq. 4.18): 

R circ = 4(1 + <z) 4 / 3 (0.500 - 0.227 LOG(q)) 4 P^, (1) 

where q is the mass ratio Md/MBH- 

Due to the viscosity, the torus at K C irc will be stretched 
and flattened into a disc on a viscous time scale (by angu- 
lar momentum transfer) . The size of the viscosity is highly 
uncertain but in the a-prescription of classical disc theory 
it is parameterized and the viscous time scale at R C irc has 
a scaling law (see FKR p. 99, eq. 5.63): 

t visc = 370(a/0.00l)- 4/5 (M BH /U) 1/4 (M d /1.5E - 8)- 3/1 ° 

x(R circ /U) 5/4 years (2) 
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where the parameters are scaled to those used in Table 
1 for the best-fit binary parameters and a — 0.001. We 
may call this the recurrence time during which a new disc 
is formed if the old one has been rapidly swallowed into 
the BH. The mass accumulated in the torus around R C i rc 
during this time equals to M accMm = t V i SC x M d - 

Surprisingly, M accum is roughly equal to the mass of 
a classical viscous disc (using a-prescription, gas pressure 
and Kramer's opacity) if the outer radius is set equal to 
R circ and 2xl0~ 7 MQ/year is used for the disc accretion. 
This accretion rate can be derived from the RXTE ob- 
servations over the past six years. The ASM light curve 
gives a time-averaged mean value of 58 counts/s (0.77 in 
the Crab-units) between 2-13 keV which corresponds to 
a total intrinsic luminosity L = 1.2 x 10 39 (d/12.5 kpc) 2 
erg/s using PCA+HEXTE fits by Vilhu et al. (2001). The 
distance d is scaled to the mean value 12.5 kpc given by 
Chaty et al. (1996) with ±1.5 kpc uncertainty. This lumi- 
nosity is slightly below the Eddington luminosity of a 14 
Mq star and corresponds to a mass accretion rate 

M obs = 2.0 x 10- 7 (r]/0.1)- 1 (d/12.5kpc) 2 MQ/year, (3) 

if n is the efficiency of converting accretion to radiation 
(L = rjM oos c 2 ). For a non-rotating black hole the max- 
imum available gravitational potential energy (per rest 
mass energy) at the innermost stable orbit is 0.06 - 0.1, 
while for an extreme Kerr-hole the efficiency may be as 
high as 0.4 (sec FKR, p. 191). 

The observed high accretion rate eats the mass from 
the torus on a timescale t ac tive = M accMm /M b s . We call 
this the 'activity time' and it is one order of magnitude 
shorter than the recurrence time = t V i SC : 

^active /^recurrence — (0.05— 0.1)(jj/0.1). 

Together these two timescales form a duty-cycle and 
their estimates are presented in Table 2. 

4. Discussion and Conclusions 

We have estimated the He-core mass (around 0.28 Mq) 
of the donor of GRS 1915+105 for the binary parameters 
given by Greiner et al. (2001b). The evolutionary expan- 
sion of the donor leads to a conservative mass transfer 
rate M d = (1.5 ±0.5) x 10~ 8 M /y which is ten times 
smaller than the accretion rate derived from the mean 
ASM light curve over the past 6 years, and using an effi- 
ciency of 0.1 to convert the mass infall into radiation: M b s 
= 2.0x(r ? /0.1)- 1 10~ 7 M /y, for a distance of 12.5 kpc. 
We propose that these two numbers determine the duty 
cycle where the active phase (as observed at present) is 
ten times shorter than the quiescent one. 

We identify the duration of the quiescent phase 
(recurrence time) as the viscous timescale at the 
circularization radius and estimate its value to be 
370x(a/0.001)~ 4 / 5 years (a = the viscosity parameter in 
the a-prescription). The corresponding active phase lasts 
28x(^/0.1)(a/0.001) -4 / 5 years and is comparable to the 



present activity phase which has already lasted for ten 
years, if the small a = 0.001 used can be justified at the 
outer edge of the disc. 

The a-parameter is highly uncertain and consequently 
so are the timescales derived. However, a comparison can 
be made with the recurrent X-ray Nova and Soft X-ray 
Transient A0620-003, using its parameters during quies- 
cence: M B h = 10 M , M d = 0.7 Mq and P = 7.75 hours 
(Tanaka and Lcwin, 1995). The accretion rate at the outer 
disc during quiescence, as derived from optical observa- 
tions (McClintock et al. 1995), equals M d = 10" 10 Mq / y 
which we identify as the mass transfer rate from the donor. 
We note that the observed X-ray luminosity during quies- 
cence implies much smaller accretion in the inner regions 
of the disc (Narayan et al. 1996). During the maximum 
outburst the accretion onto the primary black hole proba- 
bly approached the Eddington rate M b s = 10~ 7 Mg/y for 
n = 0.1, with an c- folding time of one month (Tanaka and 
Lewin, 1995). The circularization radius is 0.60 Rq lead- 
ing to t visc = 30 years (using a = 0.001 in eq. 2) which is 
briefly consistent with the two observed outbursts (1917 
and 1975) supporting the idea that the disc filling time 
is equal to the viscosity timescale at R C i rc with small a. 
The mass accumulated is small (3 x 10~ 9 Mq) and conse- 
quently the predicted active phase of A0620-003 is short 
(10 days) but of the same order of magnitude than the 
observed e-folding time. 

At present there are more sophisticated disc-models 
available than the simple a-prescription used. These in- 
clude e.g. irradiated disc-models (see King, 2000, and ref- 
erences therein). Their usage would affect the viscosity 
timescale for a fixed a but probably less than the uncer- 
tainty in a itself. 

If the mass of the donor of GRS 1915+105 is higher 
(2.4 Mq, instead of 1.2 Mq as used in the above esti- 
mate), the mass transfer from the donor will be increased 
to 5xl0~ 8 Mg/y. Further, if at the same time we in- 
crease the efficiency of the BH conversion to radiation to 
0.4 (instead of 0.1), like in the case of an extreme pro- 
grade Kerr-hole, then M d and M b s become equal. In this 
most extreme case the reasons behind the quiescent/active 
states must be searched elsewhere, e.g. in strong advection 
(ADAF) during the quiescence. We also checked that the 
hydrogen ionization zone (at around Rq) is always inside 
the circularization radius and may thus be the trigger for 
the limit-cycle instability lasting for the whole activity 
phase. 

The models in Table 1 (including uncertainties in the 
donor mass) predict bolomctric luminosities L = (50 - 
100)Lq, surface effective temperatures T e // = 3800 - 
4000 K and gravities log(g) = 1.7 - 1.9. These correspond 

to absolute K-magnitudes between —2.7 2.2 which are 

inside the limits (—2 3) given by Greiner et al. (2001a), 

but a more accurate value could properly fix the donor 
mass, and consequently its mass transfer rate. 

Another complication which should be studied is the 
possible effect of X-ray heating of the donor. Hard photons 
above 10 keV can penetrate through its photosphere into 
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the convective zone affecting its structure (Podsiadlowski, 
1991;Vilhu, Ergma and Fedorova, 1994). The mean lumi- 
nosity of GRS 1915+105 above 10 keV is roughly 2x 10 38 
erg/s (Vilhu et al. 2001) of which 0.5 - 1 per cent is cap- 
tured by the donor, assuming no screening of the disc. 
If the activity phase lasts 1/10 of the whole cycle then 
(10-30)Lq can be deposited in deep layers of the donor, 
averaged over the longer thermal timescale of the donor 
envelope, leading probably to an overestimate of the He- 
core mass and mass transfer rate. 
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Table 1. Computed parameters for the evolving donor of GRS 1915+105 using M B h = 14 M Q , M d = 1.2 M Q , P = 33.5 day 
and assuming that the donor fills its Roche lobe (first line), using the analytic methods by Webbink, Rappaport and Savonije 
(1983) for Z = 0.02. The second line gives the ranges if the donor mass is varied between (0.6 - 2.4) Mq. The He-core mass 
(Mffe), luminosity (L), radius (R), binary separation (a) and the circularization radius (R C irc) are given in solar units. The 
mass transfer from the donor (Md) is in units of solar masses per year. The predicted spectral types and absolute K-magnitudes 
(Mit) are estimated from T e ff and L using bolometric corrections and colours from Cox (2000). 



Sp M K M He L R M d a R circ 

K6 -2.6 0.28 77 21 1.5 10~ 8 108 14 

K5-M1 -2.2 - -2.7 0.26 - 0.29 50 - 100 17 - 27 5 10~ 9 - 5 10~ 8 95 - 115 12 - 18 



Table 2. Viscosity time scale t„; sc (recurrence time) at the circularization radius, the mass M accllm accumulated from the 
donor during t„; sc and the time scale t ac tive during which the BH swallows M accum with the observed mean accretion rate 
2x 10 — 7 (r;/0.1) — 1 (d/12.5 kpc) 2 MQ/year where r\ is the accretion to radiation conversion factor and d the source distance. All 
the values should be multiplied by (a/0.001)~ 4 ^ 3 where a is the viscosity parameter in the a-prescription of classical discs. The 
models are the same as in Table 1. 



t„i sc years M accum /M sun t ac tive / (f?/0.1) years 

370 5.5 x 10~ 6 28 

200 - 700 (3.7 - 9.0) x 10~ 6 20 - 45 



